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In a data sample of approximately 1.3 fb _1 collected with the DO detector between 2002 and 2006, 
the orbitally excited charm state 75^(2536) has been observed with a measured mass of 2535.7 ± 
0.6 (stat) ± 0.5 (syst) MeV/c 2 via the decay mode Bg — > Z?7i(2536) /i + v^X . A first measurement is 
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made of the branching ratio product Br(b — > D7i(2536)/i + v ll X)-Br{D~ 1 — > D*~ K%). Assuming that 
D ~l (2536) production in semileptonic decay is entirely from B®, an extraction of the semileptonic 
branching ratio Br(B® — > Djj (2536)/x + ^ M X) is made. 

PACS numbers: 13.25.Hw,14.40.Lb 



Semileptonic B° s decays into orbitally excited P-wave 
strange-charm mesons (D**) are expected to make up 
a significant fraction of J5° semileptonic decays and are 
therefore important when comparing inclusive and exclu- 
sive decay rates, extracting CKM matrix elements, and 
using semileptonic decays in B® mixing analyses. For 
B meson semileptonic decays to heavier excited charm 
states, more of the available phase space is near zero 
recoil, increasing the importance of corrections in heavy- 
quark effective theory (HQET) [l[ , effectively tested here. 

D** mesons (also denoted D sJ ) are composed of a 
heavy charm quark and a lighter strange quark in an 
L = 1 state of orbital momentum. In the heavy-quark 
limit, the spin sq of the heavy quark and the total angu- 
lar momentum, j q = s q + L of the light degrees of free- 
dom (quark and gluons), are separately conserved and 
the latter has possible values of j q — \ or | . The sur- 
prisingly light masses of the j g = \ states: D* (2317) and 
Z3 s i(2460) @], plus the observation of new D s j states 
deepens the need for a better understanding of these D** 
systems since they may be quark molecular states, a new 
and very different arrangement of quarks. 

In our decay of interest, the j q = | angular momentum 
can combine with the heavy quark spin to form the J p = 
1 + {D s i) state which must decay through a D-wave to 
conserve j q — |. The £^(2536) is expected to decay 
dominantly into a D* and K meson to conserve angular 
momentum. 

In this Letter we present the first measurement of 
semileptonic B® decay into the narrow Df x (2536) state. 
This state is just above the D* K s mass threshold and 
has been observed previously [3| • Events compatible with 
the decay chain b_ -> D~ L (2536)/i+^X, ^(2536) -» 
D*~K%\ D*~ -> £>°7r-, K% -> 7T+7T-, D° -► K+tt~ are 
reconstructed. Charge conjugate modes and reactions 
are always implied in this Letter. 

Assuming that Z?7i(2536) production in a semilep- 
tonic decay is entirely from B®, the branching ratio 
Br{B° s — » Z)^(2536)/i + z/ Al X) can be determined by nor- 
malizing to the known value of the branching fraction 
Br(b -> D*-fi+v^X) = (2.75 ± 0.19)% @ to avoid un- 
certainties in the &-quark production rate. This semilep- 
tonic branching ratio includes any decay channel or se- 
quence of channels resulting in a D* and a lepton (muon 
in our case), and all b hadrons, and therefore includes 
the relative production of each b hadron species start- 
ing from a b quark. Since the final state of interest, 
D~ L (2536) — > D*~K S , is reconstructed from a D* and 
a K s , the selection is broken up into two sections: one 
to reconstruct the D* with an associated muon, coming 



dominantly from B meson decays resulting in a number 
of candidates, iVo*^, and then the addition and subse- 
quent formation of a vertex of a Kg with the D* and 
muon, resulting in Njj 3l candidates. To find the branch- 
ing ratio, the following formula is used: 

f(b^B° s ).Br(B° s ^D^+^Xy 

■Br(D 7l - D*-K° s ) = Br(b - D^^X) ■ 

e(b^D*fi) J_ 
e(B° s -> D slf i -» D*fi) e K o- { > 

The input f(b—>Bg) [5[ is the fraction of decays where 
a b quark will hadronize to a B® hadron. e K o is the 
efficiency in the signal decay channel to reconstruct and 
make a vertex with a K s to form a Z? s i(2536), given that 
a D* and a muon have already been reconstructed. Later 
we will identify the ratio of efficiencies as i?fjf = £(B® — ► 
D sl ^ -> D*fx)/e(b -> D*n). 

The DO detector [6j and following analysis [7] are de- 
scribed in more detail elsewhere. The main elements rel- 
evant to this analysis are the silicon microstrip tracker 
(SMT), central fiber tracker (CFT), and muon detector 
systems. 

This measurement uses a large data sample, corre- 
sponding to approximately 1.3 fb _1 of integrated lumi- 
nosity collected by the DO detector between April 2002 
and March 2006. Events were reconstructed using the 
standard DO software suite. To avoid lifetime biases 
compared to the MC simulation, the small fraction of 
events were removed that entered the sample only via 
triggers that included requirements on impact parame- 
ters of tracks. 

To evaluate signal mass resolution and efficiencies, 
Monte Carlo (MC) simulated samples were generated for 
signal and background. The standard DO simulation and 
event reconstruction chain was used. Events were gen- 
erated with the PYTHIA generator Q and decay chains 
of heavy hadrons were simulated with the evtgen de- 
cay package Q. The detector response was modeled by 
GEANT [lfj . Two background MC samples were also gen- 
erated: a cc sample, and an inclusive 6-quark sample 
containing all b hadron species with forced semileptonic 
decays to a muon. In both cases, all events containing 
both a D* and a muon were retained. 

B mesons were first selected using their semileptonic 
decays, B — > D*~ /i + X. At this point in the selection, the 
D*+n sample is dominated by B® — » D*~ [i+v^X decays. 
For this analysis, muons were required to have hits in 
more than one muon layer, to have an associated track 
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in the central tracking system, and to have transverse 
momentum pj. > 2 GeV/c, pseudorapidity < 2, and 
total momentum p^ > 3 GeV/c. Two oppositely charged 
tracks with p T > 0.7 GeV/c and |r/| < 2 were required 
to form a common D° vertex which were then combined 
with a muon candidate to form a common decay point 
following the procedure described in Ref. For each 
_D°/i + candidate, an additional soft pion was searched 
for with charge opposite to the charge of the muon and 
Pt > 0.18 GeV/c. The K~ and tt + from the decay of the 
D° were both required to have more than five CFT hits. 
To reduce the contribution from prompt cc production, 
a requirement was made on the transverse decay length, 
L xy , significance of the D*/j, vertex of L xy /o~(L xy ) > 1. 
After these cuts, the total number of D* candidates in 
the mass difference, M{D*) - M(D°), peak of Fig. ED is 
N D ^ = 87506 ±496 (stat). 
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FIG. 2: Invariant mass of D*Kg with an associated muon. 
Shown is the result of the fit of the D* K% mass with the 
function described in the text. 
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FIG. 1: The mass difference M(D*) - M(D°) for events with 
1.8 < M(D°) < 1.95 GeV/c 2 and an associated muon. The 
number Nd'^ was defined as the number of signal events in 
the mass difference range of 0.142-0.149 GeV/c 2 . 



D sl (2536) candidates were formed by combining a 
D* candidate with a Kg. D* candidates were first se- 
lected by requiring the mass difference M(D*) — M (D°) 
to be in the range 0.142-0.149 GeV/c 2 . The two tracks 
from the decay of the K s were required to have opposite 
charge and to have more than five hits in the CFT detec- 
tor. The pt of the K s was required to be greater than 
1 GeV/c to reduce the contribution of background K s 
mesons from fragmentation. A vertex was then formed 
using the reconstructed K s and the D* candidate of the 
event. The decay length of the K s was required to be 
greater than 0.5 cm. To compute the D^(2536) invari- 
ant mass, a mass constraint was applied using the known 
D*^ mass 0| instead of the measured invariant mass of 
the i^7r7r system. Finally, the invariant mass of the re- 
constructed -D^(2536) and muon was required to be less 
than the mass of the B° s meson [j| . 

The signal model employed for the fit to the D*Kg in- 
variant mass spectrum was a relativistic Breit-Wigner 



convoluted with a Gaussian function, with the reso- 
nance width fixed to the value 1.03 ± 0.05 (stat) ± 
0.12 (syst) MeV/c 2 measured by the BaBar Collabo- 
ration [12j and a Gaussian width determined to be 
2.8 MeV/c 2 from MC simulation of the signal. The MC 
width value was scaled up by a factor of 1.10 ± 0.10 to 
account for differences between data and MC resolution 
estimates. The unbinned likelihood fit used an exponen- 
tial function plus a first-order polynomial to model the 
background with a threshold cutoff of M(D*) + M{K° S ). 
The fit, shown in Fig. O gives a central value for the 
mass peak of 2535.7 ± 0.7 (stat) MeV/c 2 , a yield of 
jVd s i = 45.9 ± 9.1 (stat) events, and a significance of 
6.1cr for the background to fluctuate up to or above the 
observed number of signal events. 

The efficiencies used in Eq. [1] are estimated using the 
MC simulation, after implementing suitable correction 
factors to ensure proper modeling of the underlying b- 
hadron pt spectrum, as well as trigger effects. An event- 
by-event weight, applied as a function of the generated 
Pt of the B s , was determined by comparing the gener- 
ated pt(B) in MC with the pt distribution of fully re- 
constructed B + — > J/tpK + candidates in data collected 
primarily with a dimuon trigger |13| . Most events for 
this analysis were recorded using single muon triggers, 
and an additional weight was applied as a function of 
Pt(m) to further improve the simulation of trigger effects. 
Reweighted MC events were used in the determination of 
efficiencies described below, and indicated uncertainties 
are due to MC statistics. 

Using the MC sample of inclusive b — ► D*fj,X events, 
specific major decay modes were identified. Efficien- 
cies for each of these decay modes to pass the D*fi se- 
lection, including the efficiency to reconstruct the soft 
pion from the D* , were then determined. The pre- 
dicted fraction Fj of each channel contributing to the 
D* p sample before further cuts was found following a 
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procedure similar to that given in Ref. 14]. The ef- 
ficiency ej for each channel was found and a weighted 
sum was calculated, giving an estimated total efficiency 
for reconstruction of e(6 -► D*p) = (5.88 ± 0.80)%, 
where the uncertainty is dominated by the MC statistics 
used to find e*, and uncertainties on external inputs [B| 
used to estimate Fi. Applying the same cuts for recon- 
structing the D*n for the signal channel, the efficiency 
e(B° s -> D sl [i -> D*p) = (3.20 ± 0.02)%, results in a 
ratio of efficiencies of i?f,°, n = 0.547 ± 0.075. 

The signal MC sample was used to determine the ef- 
ficiency to reconstruct 1571(2536) — > D*~K S given a re- 
constructed D* fx as a starting point. This efficiency is 
hence effectively that of reconstructing a Kg — » 7r + 7r~ 
and forming a vertex with the D*u, and includes the 
branching ratio Br(Kg — > Tr + n~) [5j for ease of use in 
calculating the branching ratio product. The reconstruc- 
tion efficiency was found to be e^o = (10.3±0.4)% where 
the uncertainty is due to MC statistics. 

The process cc — > D*~ [i + v ^X can contribute to No* ^ 
since a D* meson can come from the hadronization of the 
c quark, and the muon can come from the semileptonic 
decay of the hadron containing the c quark. To determine 
the number of events in our signal reconstructed from a 
prompt D* , a comparison was made of the decay length 
significance distribution observed in the data with the 
same distribution predicted by MC for b — > D* [iX and 
any excess at shorter significances was interpreted as cc 
contribution. For the decay length significance cut used 
in the analysis, L xy /a(L xy ) > 1, the fraction of No-^ 
from cc production was estimated to be (3.9 ± 2.5)%. A 
check using a prompt cc MC sample results in a consistent 
estimate. The value of Nd*h was corrected downward 
accordingly. 

The contribution from cc production to No 3l where 
one charm quark hadronizes directly to a -D s i(2536) and 
the other decays directly to a muon was estimated to 
be negligible using relative production ratios and spin- 
counting arguments 
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Systematic uncertainties for the branching ratio prod- 
uct are summarized in Table U and discussed below. 
The uncertainty in the normalizing branching ratio [|| 
Br(b — > D*fiX) was taken as a systematic uncertainty. 
For determining Np*^, the signal and background model 
parameters were varied in a correlated fashion and a 
systematic uncertainty was assigned. The estimated cc 
production contribution was varied by the indicated un- 
certainty. In the determination of Nd s1 , the functional 
forms of the signal and background models were varied in 
a number of ways to determine the sensitivity of the can- 
didate yield. In addition, the scaling of the widths was 
varied by ±10% to check the sensitivity to uncertainty in 
mass resolution. 

By comparing the £>t(m) distribution for the signal us- 
ing the default ISGW2 decay model to the HQET 
semileptonic decay model [9J, a weighting factor was 



found and applied to the fully simulated signal MC 
events, and the efficiency determined again. The dif- 
ference observed was assigned as a contribution to the 
systematic uncertainty of e^-o and Rp™. 

When estimating e K o , the uncertainty due to model- 
ing of the b hadron px spectrum was derived by using 
an alternate weighting technique. The cuts on the px 
and decay length of the K s were varied and a system- 
atic uncertainty on the efficiency due to this source was 
also assigned. Discrepancies in track reconstruction ef- 
ficiencies between data and MC in low-pr tracks were 
accounted for by assigning a systematic uncertaint y t o 
each of the pion tracks in the K s reconstruction IT], ll8( . 

The uncertainty in R%, is due to a combination of 
MC statistics and uncertainties in PDG branching ratio 
values and production fractions, f(b — > b hadron). The 
uncorrelated systematic uncertainty is given in Table [U 

The estimated systematic uncertainties were added in 
quadrature to obtain a total estimated systematic un- 
certainty on the branching ratio product of 16.8%. The 
branching ratio product was determined to be: 



Br(B° s - D-^+^X) ■ Br(D 7l - D* 

= [2.66 ± 0.52 (stat) ± 0.45 (syst)] x 10~ 4 . 



TABLE I: Estimated systematic uncertainties. 



Source 



Systematic uncertainty 



Br(b - 

Nd s1 

risen 
-fin* 



D*fiX) 



6.9% 
2.9% 
5.5% 
11.0% 
8.6% 



Total 



16.8% 



To assess the systematic uncertainty on the mass mea- 
surement, the same variations of the Z? s i(2536) mass sig- 
nal model, as well as background functional form, were 
applied as described above. The mass values used for 
the mass constraints on the decay products were varied 
within their PDG uncertainties and were also set to the 
DO central fit values. Ensemble tests indicated that the 
statistical error is correct. From the observed variations, 
a total systematic mass uncertainty of 0.5 MeV/c 2 was 
taken, for a mass measurement of: 

m(D sl ) = 2535.7 ± 0.6 (stat) ± 0.5 (syst) MeV/c 2 . 

This measured mass value is in good agreement with the 
PDG average value of 2535.34 ± 0.31 MeV/c 2 

To allow comparison of this measurement to theoreti- 
cal predictions, the semileptonic branching ratio alone as 
shown in Table Hi] is extracted by taking the hadroniza- 
tion fraction into B° as f(b -> B° s ) = 0.103 ± 0.014 @ 
and also assuming that i?r(D s i(2536) — > D*Kg) = 
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0.25 [9(. This is the first experimental measurement of 
this semileptonic branching ratio and is compared to a 
number of theoretical predictions 0, [2(| of the ex- 
clusive rate in Table [IT] The systematic uncertainty on 
this quantity is as described earlier, and the error la- 
beled "(prod, frac.)" is due to the current uncertainty 
on f(b — > B®). The first two theoretical predictions in- 
clude relativistic and corrections, while the third 
does not. The result is found to be consistent within un- 
certainties with the first two theoretical predictions, and 
demonstrates the need for such corrections. 

TABLE II: Experimental measurement compared with vari- 
ous theoretical predictions. 

Source Br(B" — > DJ 1 (2536)fi + u^XY 

This result [1.03 ± 0.20 (stat) ± 0.17 (syst) 

±0 14 (prod.fr ac.)]% 

Theoretical Predictions Br(B" -> #'(2536)^+^) 

ISGW2 [1] (0.53 ± 0.27)% 

Relativistic Quark Model & 

1/rng corrections [19] (1.06 ±0.16)% 
Non-rel. HQET and ISGW [20] 0.195% 



In summary, using 1.3 fb _1 of integrated luminosity 
collected with the DO detector, a first measurement of the 
semileptonic B® decay into the narrow 15^(2536) state 
has been made and compared with theory. In addition, 
the mass of the £>^ (2536) was measured and found to be 
in good agreement with the PDG value. 
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